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ABSTRACT 



The ensemble optical/ultraviolet variability of narrow-line Seyfert 1 (NLS1) 
type active galactic nuclei (AGNs) is investigated, based on a sample selected 
from the Sloan Digital Sky Survey (SDSS) Stripe-82 region with multi-epoch 
photometric scanning data. As a comparison a control sample of broad-line 
Seyfert 1 (BLS1) type AGNs is also incorporated. To quantify properly the in- 
trinsic variation amplitudes and their uncertainties, a novel method of parametric 
maximum-likelihood is introduced, that has, as we argued, certain virtues over 
previously used methods. The majority of NLSl-type AGNs exhibit significant 
variability on timescales from about ten days to a few years with, however, on 
average smaller amplitudes compared to BLSl-type AGNs. About 20 NLSl- 
type AGNs showing relatively large variations are presented, that may deserve 
future monitoring observations, for instance, reverberation mapping. The av- 
eraged structure functions of variability, constructed using the same maximum- 
likelihood method, show remarkable similarity in shape for the two types of 
AGNs on timescales longer than about 10 days, which can be approximated by a 
power-law or an exponential function. This, along with other similar properties, 
such as the wavelength-dependent variability, are indicative of a common dom- 
inant mechanism responsible for the long-term optical/UV variability of both 
NLS1- and BLSl-type AGNs. Towards the short timescales, however, there is 
tentative evidence that the structure function of NLSl-type AGNs continues de- 
clining, whereas that of BLSl-type AGNs flattens with some residual variability 
on timescales of days. If this can be confirmed, it may suggest that an alter- 
native mechanism, such as X-ray reprocessing, starts to become dominating in 
BLSl-type AGNs, but not in NLS1-, on such timescales. 
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Introduction 



Variations of optical/ultraviolet (UV) luminosity on timescales from weeks to years 
are characteristics of active galactic nuclei (AGNs), and the study of this phenomenon is 
a powerful tool for constraining models of black hole accretion. Extensive observational 
investigations in the past, albeit with mostly small samples, have revealed dependence 



of the variability amplitude o n various observec 



luminosity, and redshift (e.g. 



di Clemente et al. 



quantities such as wavelength, time-lag, 



1996 



de Vries et al. 



2005 



2011 



Welsh et al 



Meusinger et al. 



20111 ). that was confir med by using vast quasar s from the Sloan Digital 



Sky Survey (SDSS) with repeated scans fjVanden Berk et al. 



200 1 



Zuoetal 



2012J). With 



the advent of black hole mass estimation based on single-epoch spectroscopic data on AGN 



new c orrelations have b een suggested of the variations with black hole mass (IWold et al. 



2007 



Bauer et al. 



20091 ). and perhaps more fundamen tally, with the 



claimed independent ly by a few groups very recently (IWilhite et al. 



Edd i ngton ratio as 



20081: 



Ai et al. 



2010 



MacLeod et al. 



20101 ) . The possible dependence of variability on the fundamental physical 
parameters is remarkable, as it may provide new insights into the variability mechanisms. 
For instan ce, such a dependence favors models of disc instability o r variations in accretion 



rate (e.g., 



Kawaguchi et al. 



1998 



Li fc Ga 



o 



2008 



such as Poisson process (e.g. 
gravitational microlensing (e.g. 



Terlevich et al. 



Hawkins 



1992 



Liu et al. 



2008 ) over other mod els 



Torricelli-Ciamponi et al. 



20001) and 



19931 ). However, the statistical significance of 



these results is not yet sufficiently high, and further confirmation is still needed. 

While considerable progresses have been made in the past decades, previous variability 
studies were focused mostly on AGNs of the broad-line Seyfert 1 (BLS1) type, including 
quasars. Narrow-line Seyfert 1 (NLS1) type AGNs, as a subclass of broad emission-line 
AGNs (BLAGNs) with extreme properties, have received little attention, however. This 
may partly be ascribed to the small number of such objects known in the past. NLSl-type 
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AGNs are commonly defined a s having small widths of t 



(FWHM(if/3) < 2000 km s" 1 , 



Osterbrock fe Pogge 



le optical permitted emission-lines 



1985|). In addition, NLSl-type AGNs 



see Komoss a 2008 for a review) often show strong optical Fe II e mission and weak 



[O III] lines flGaskell 



spectral slop e 



1999 



Grupe 



1985 



Boiler et al 



Goodrich 



1996 



1989 



Veron-Cettv et al. 



Wang et al 



20011 ) . a steep soft X-ray 



19961 ). and strong X-ray variability ( iLeighly 



20041 ). Such extreme properties make them cluster at o ne extreme end 



of th e strongest AGNs correlati on space, namely eigenvector 1 (EV1, 



1992 



Sulentic et al. 



2000 



Boroson fe Green 



20071 ). which is believed to be driven by the Eddington ratio 



EddJ 



NLSl-type AGNs are commonly thought to have small black hole masses and thus 



high accretion rates (at a substantial frac tion of or close to the 



argued extensively in the 



Collin fe Kawaguchi 



2001 



iteratures (e .g. 



Grupe 



Mineshige et al. 



2000 



Eddington rate), as 



Sulentic et al. 



2000 



20041 ) . As such, their accretion is suggested to likely 
proceed via the so-called slim disk where a significant fraction of viscosity-generated photons 
are trappe d and advected into the black hole, rather than the standard Shakura-Sunyaev 



thin disk (Abramowicz et al. 



1988 



Wang fe Netzer 



20031 ). As the accretion disk is 



responsible for the UV/optical radiation in AGNs, different variability patterns might be 
anticipated for NLS1-, compared to normal BLSl-type AGNs. 

There have been a few previous variability studies on NLSl-type AGNs, all of which 



were monitoring programs on relati vely short timesca 



individual objects or small samples (jYoung et al 



1999 



es (e.g. intra-nigh t or days) for 



Miller et al. 



20001 ). The only one 



NLS1 with long-term (m ore than ten years) extensive photometric monitoring is Ark 564 



( iDoroshenko et al. 



20061 ). as far as we know, and it was found that the NLS1 has a low 



variability amplitude. Based on monitoring results of a small sample, Klimek et al. (2004) 
first suggested that NLSls, as a class, are less variable than non-NLSls on long timescales 
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(weeks to years). This is consistent with the later suggested anti-correlation with Eddington 
ratio, if NLSl-type AGNs have indeed high accretion rates. 

In fact, the low black hole mass scenario of NLSl-type AGNs is yet subject to active 
debates (see komossa 2008 for a review), albeit some favoring independent pieces of indirec 



evide nce, such as rapid X-ray variability ( iGierliriski et al. 



2008 



Ponti et al. 



2010 



Ai et al. 



201ll ). One important question concerned is whether NLS1 AGNs conform the same relation 
of the size of the broad line region (BLR) and luminosity as found in normal Seyferts 
and quasars (Kaspi et al. 2000). The question arises as NLS1 AG Ns with reverberatio n 



mapping (RM) o bservations are rare, only limited to a few cases ( [Peterson et al 



Kaspi et al. 



2004 ; 



20051 ). This situation may partly be ascribed to the lack of known NLSl-type 
AGNs with large variability. A search of more such objects would be essential for future 
RM observations of NLS1 AGNs. 

We have carried out a systematic study on the optical/UV variability of AGNs and 
quasars using a large sample with multi-epoch scanning data available from the SDSS 
Stripe-82, including both BLS1- and NLSl-type AGNs. As the firs t results, variou s 



correlations with AGN variability have been tested and presented in 



Ai et al. 



(120101 ). It 



was found that the variability amplitude is significantly correlated with the A GN EV1 



param eters, particularly with Lboi/^Edd as the strongest, albeit large scatters (lAi et al. 
2010|). This is also manifested by the observed averagely smaller variations of NLS1- (having 
higher Lb i/-^Edd) compared to the BLSl-type AGNs. However, these results are based on 
time-averaged variation amplitudes, and confirmation by structure function analysis, that 
is a more robust technique, is needed; moreover, the time-lag dependence can be explored, 
that may shed light on the variability difference between the two types of AGNs. 

This is the second paper of our variability study based on the same SDSS Stripe-82 
samples, which focuses on the ensemble variability properties of NLSl-type AGNs. The 
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BLSl-type subsample is also analyzed in the same way as a control sample for comparison. 
These mainly concern the demography of variable objects, wavelength dependence, and 
variation of the spectral shapes, as well as the dependence on time-lag (structure function). 
We also improve the measurement of variation amplitudes up on those in the first paper (and 
in the literature) by introducing the Maximum-likelihood method. In addition, we present 
a full account of the data analysis, mainly photometric calibration and reliability check, for 
which only a concise description was given in the first paper given the limited length of a 
letter. The paper is organized as follows. In Section 2 we describe the sample selection and 
photometric data analysis, and in Section 3 the quantification of variability. The results 
are presented in Section 4, and their implications are discussed in Section 5, followed by a 
summary in Section 6. We use a A-dominated cosmology with Hq = 70 km s" 1 Mpc^ 1 , 
Q m = 0.3, and Qa = 0.7 throughout the paper. 



2. Sample and Data 



2.1. NLS1- and BLSl-type AGN samples 



Our NLSl-type AGNs are actually taken from the sample of 



Zhou et al. 



(120061 ). which 



was selected from the SDSS Data Release 3 catalogs. The procedures of the optical spectral 



analysis and me asurement o 



previous papers (IZhou et al 



emi s sion line paramet ers were described in detail in our 



2006 



Dong et al. 



20081 ). Only objects with redshifts z < 0.8 



were included to cover the H/3 and [O III] emission lines, with the broad component of Ha 
or H/3 detected at > 10 a confidence level. We further select objects for variability study 
using the following criteria: 



located within the SDSS Stripe-82 region (R.A. > 310° or < 59° and -1°.25 > deck < 
1°.25); 



classified as 'STAR' in all of the five bands by the SDSS photometric pipeline^ to 
minimize possible contamination from host galaxy starlight; 



radio-quiet or not detected in the FIRST radio s urvey 



(Faint Images of the Radio Sky at 20 cm survey, 



, as constrain ed from the FIRST 



Becker et al. 



19951 ) radio survey. This 



is to eliminate possible 'contamination' from blazar-type jet emission, w 



rich origina tes 



from a distinct (non-thermal) radiation process and is o: 



2003 . 



2007 



Yuan et al. 



2008 



Abdo et al. 



2009 



Liu et al. 



ten variable 



2010 



Zhou et al. 



Jiang et al. 



20121 ). 



The above selections result in 58 NLSl-type AGNs. For a comparative study we also 
select a BLSl-type AGN sample in the Stripe-82 region, i.e., with the broad Ha or H/3 
linewidth (full width at half maximum, FWHM) greater than 2,200 km s _1 , the dividing 
line used to separate NLS1- and BLSl-type AGNs following our previous work (Zhou et 
al. 2006, see also Gelbrod et al. 2009). The sample selections and spectral a nalysis are 



perfor med in the same way as for NLSl-type AGNs above, and is described in 



Dong et al. 



( 120081 ). The same criteria as above are also applied to select star-like, radio-quiet objects. 
Furthermore, To ensure that the two samples match each other on the redshift-luminosity 
distribution, a subsample of the resulting BLSl-type AGNs is randomly extracted in a 
way that they mimic the z — Mi distribution of the NLSl-type sample, where Mi is i-band 



x The SDSS photometric pipeline differentiates resolved ('GALAXY') and unresolved 
('STAR') sources based on the difference between PSF- and model-magnitudes, and it was 
found that a simple cutoff at 0.145 mag works at the 95% confidence level for sources with 



r < 21 mag (IStoughton et al. 



20021 ). 



2 Defined as -R1.4 < 1, where the radio-loudness -R1.4 is defined as the logarithm of the flux 
density ratio of radio emission at 1.4 GHz to the optical emission at the SDSS g-band. 
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absolute magnitude In this process we try to retain the size of the NLSl-type sample 
(with only a few outliers discarded) and prune the BLSl-type sample with a BLSl-type 
to NLSl-type ratio consistent with 2:1. This results in two final working samples of 55 
NLSl-type and 108 BLSl-type AGNs, which are statistically compatible with each other o n 



the r edshift-luminosity plane (Figured]). The 2-D Kolmogorov-Smirnov test (IPress et al. 
19921 ) yields a chance probability of 0.43 that the two samples have the same z — M\ 
distribution. 

The two working samples can be considered as optically and ho mogeneously select ed 



with reliab 



Dong et al. 



measured continuum and emission line parameters (see 



Zhou et al 



2006: 



20081 ). They consist of mostly quasars with Mj < —23 mag and z ~ 0.3-0.8. 
All of the NLSl-type objects meet the conventional extra [O III]/H/3< 3 criterion for 
NLSls (see zhou et al. 2006, for details). They also show, compared to the BLS1 sample, 
other properties characteristic of NLSls, e.g., stronger Fe II emission, smaller black hole 
masses and higher Eddington ratios (Figure[2]). The strength of the Fe II emission lines 
is measured by the intensity ratios to H/3, i.e., i?457o=Fe II(AA4434 — 4684)/H/3 b (where 
Fe II (AA4434 - 4684) denotes the Fe II multiplets flux integrated over 4434-4684 A, and 
H/3 b the flux of the broad component of H/3, see Zhou et al. 2 006 for details ). The black 



Greene fc Hoi fl2005h from 



hole masses, Mbh, are estimated using the formalism given by 
the br oad H/3 FWHM and the 5100 A luminosity (L5100), which are taken from 



Zhou et al 



(120061 ). or measured in the same way for the BLS1 sample. To calculate the 



ratio , L bo i/L Edd , the bolometric luminosities are estimated as L bol = 9L 510 o (jElvis et al. 



Eddington 



1994J). 



3 A power law continuum fx oc X ax with a slope of ot\ = —1.5 is assumed in calculating the 
i-band abs olute magnitude. Th e observed magnitudes have been corrected for the Galactic 



extinction (ISchlegel et al. 



1998|). 
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2.2. Data and photometric calibration 



The SPS S imaged the 



u, g, r, i and z (jGunn et al. 



sky i n five broad photom etric bands simultaneously, namely, 



2006 



Fukugita et al 



19961 ). The integration time is 54.1s in each 



band and the limiting magnitude reaches ~ 23 mag in the r band. The phot ometric syste m 
is ba s ed on the AB-system with zero-point u ncertainty of ~0. 02-0. 03 mag (ISmith et al. 



2002 



Fukugita et al. 



1996 



Ivezic et al. 



20041 ). The astrometric accuracy is better than 0."1 



20031 ) . We use the point-spread 



for sources brighter than 20.5 mag in r band ( jPier et al. 
function (PSF) magnitudes in this work. 

The SDSS Stripe-82 survey covers the region from a = 59° to 310° and 5 = —1.25° 
to 1.25°. During the SDSS-I phase (~ 2000-2005) the region was repeatedly observed 
and the central part of the strip e has been scanne d at a cadence of typically 10-20 times 



( lAdelman-McCarthy et al. 



2007 



Sesar et al. 



20071 ). These observations were performed 



under generally photometr ic conditions and t he data were well calibrated using the 



PHOTO software pipeline fjLupton et al 



20021 ). This region was later repeatedly scanned 



over the course of three 3-month campaigns (Sep.-Nov.) in the successive three years 
in 2005-2007, known as the SDSS Supernova (SN) Survey. In this work we use the 
photometric data acquired dur ing both the SDSS-I phase from Data Release 5 (DR5, 



Adelman-McCarthy et al. 



20071 ) and the SN survey during 2005 (SN-2005). 



Observations in the SN survey were sometimes performed in non-photometric 
conditions. At the time when this work was started only the un-calibrated source catalogs 
were available. Thus possible photometric zero-point offsets in the SN survey data need 
to be determined, by calibrating against the DR5 magnitudes. To do this, we use stars in 
the same fields as 'standards', assuming that the vast majority of stars do not vary. The 
detailed procedures of re-calibration are described in Appendix [A] The resulting overall 
(systematic and statistical) photometric errors of the calibrated SN survey magnitudes 
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have a median of ~ 0.03 mag for the g, r, and i bands, and ~ 0.04 mag for the u and z 
bands, which are comparable to those of the SDSS-I DR5 data. To check the reliability 
of our photometric calibrati on, we examin e the calibrated SN survey data of 14 Landolt 



photometric standard stars (ILandolt 
show detectable variability. 



1992) locating in Stripe-82; none of them is found to 



3. Quantification of variations and their uncertainties 

We construct light curves for each object in the five SDSS bands separately using the 
DR5 and the re-calibrated SN-2005 Survey data. Objects observed in the same frames 
observed at different epochs are matched by using a matching radius of 1". To eliminate 
effects of observation condition, only frames with good image quality (flagged as good or 
acceptable) are used. For each object there are typically ~27 observations (about 14 from 
the SN-2005 survey) spanning ~5 years. It should be noted that all sample objects have 
similar sampling patterns, mostly being observed repeatedly in three months (Sep. -Nov.) 
each year during 2000-2005. 

First we use the chi-square method to test the significance of variability for 



individual objects against t 



Bevington fc Robinson 



re null hypothesis of no variation. The statistic x ( e -g-; 



19921 ) is defined as 



x2 = E (m,-W)^ (1) 

i=l ^ 

where m; is the magnitude of the i-th measurement with error £j (the overall error 
including both statistical and systematic errors), and (m) the weighted mean of a total of 
N measurements. We consider a source to be variable only if the probability level against 
the null hypothesis P < 0.1%. 

The amplitude of variability of an object is commonly measured by the width of its 



- 12 - 



magnitude distribution. However, the magnitude distribution as observed does not represent 
the intrinsic variability, but rather is broadened by the effect of measurement errors. To 
assess the intri nsic magnitude distribution, we first take th e common practice as used in the 



literature (e.g., 



Vanden Berk et al. 



2004 ; 



Sesar et al. 



20071 ). by estimating the contribution 



of measurement errors to the observed scatter (root-mean- square, rms) and then subtract 
it from the latter (See Eq. B2 in Appendix[B] for a detailed account). This estimation is 
most accurate when the magnitude uncertainties £j have the same or very similar values. 
However, this may not be the case for the SDSS SN survey data that taken under variance 
of observing conditions. Moreover, for such estimation and its variants, it is difficult to 
quantify the uncertainties of the estimated intrinsic variability amplitudes. 

Alternatively, we introduce a parametric maximum-likelihood method, in which the 
intrinsic magnitude distribution of an object is parameterized as a Gaussian with a mean 
(m) and a standard deviation a m ; the standard deviation a m is used as a measure of the 
amplitude of intrinsic variability. This method can take into account the uncertainty of 
each individual measurement (assumed to be Gaussian distributed). Moreover, it can 
also be used to quantify the confidence intervals for interesting parameters, such as cr m , 
(Figure[3]). The method had been used in a similar way in the literature to derive the 



intrinsic distribu 


;ion 


of 


Maccacaro et al. 


1988 


)• 



19881 ). The description of this method and a comparison with the one 
above are in Appendix[B] The two methods are found to be in excellent agreement. Hence 
in this paper we use the amplitudes derived from the maximum-likelihood method since 
their confidence intervals can be estimated for each object. 
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4. Results 

4.1. Demography of Variability 

Here we quantify the demography of the variability We consider an object to be 
variable if the x 2 null- variability probability P < 0.1%, which corresponds to xt — 3 for a 
degree of freedom of 27, roughly the mean number of data sampling of our objects. Tabled] 
lists the fraction of variable objects for both the NLS1- and BLSl-type AGN samples. We 
find that the vast majority of NLSl-type AGNs show variations on timescales of years 
at the significance level P < 0.1%. Most of the variable objects varied at a level, cr m 
> 0.05 mag, larger than the typical photometric errors. This is clearly demonstrated in 
FigureHl which shows the reduced chi-square versus the variability amplitude in the r band 
as example. We can see that for NLSl-type AGNs the intrinsic variability amplitude <r m 
ranges from ~ 0.25 mag down to a few per cent at xl = 3; the latter can be considered 
as the sensitivity of the variability detection limited by the current data quality. As can 
be seen, the optical/UV emission of NLSl-type AGNs are indeed variable on timescales of 
years, with mean variability amplitudes around 0.07-0.11 mag, depending on bands. 

The accumulative distributions of the variability amplitudes <r m in the five SDSS bands 
are shown in FigureO Several results can be inferred immediately. First, most of both the 
NLS1- and BLSl-type AGNs are variable on timescales of years, and there is a large scatter, 
stretching over about one order of magnitude, in the variability magnitude. Second, the 
NLS1- tend to vary with smaller amplitude compared to the BLSl-type AGNs, and this 
trend is most significant in the blue bands. This result is in line with our previously found 
correlations involving variability, i.e., the narrower line width, the stro nger Fe II emission, 



20101). Third, For both 



the higher Eddington ratio, and the smaller variability flAi et al. 
NLS1- and BLSl-type objects the variations are larger in the bluer band than the redder 
band; this is discussed in detail in Section 4.2. 
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Variable NLSl-type AGNs with relatively large amplitude are of particular interest 
and suitable for future spectroscopic monitoring and reverberation mapping observations. 
We list objects with r band <r m > 0.1 ma g in TableEfl includ ing some of the continuum and 



emission line parameters, as taken from 



Zhou et al. 



(120061 ). We note that most of these 



are typical NLSl-type AGNs with strong Fe II emission. Their light curves are shown in 
FigureEJ 



4.2. Wavelength dependence and spectra variability 

Wavelength-dependent va riability has long been found for BLAGN and QSO samples 



[e.g., 



Vanden Berk et al. 



20041 ) . Here we examine this property for NLSl-type AGNs using 
our large data set. Due to redshift, a given filter samples different wavelength ranges for 
objects with different redshift. We approximate the derived amplitude er m for a filter band 
as that at its effective wavelength Ao- For an object of redshift z, Xq corresponds to a 
wavelength A = Ao/(l + z) in the object's rest frame. Thus for each object the variability 
cr m at five rest-frame wavelengths are sampled corresponding to the five SDSS bands, we 
divide the overall rest wavelength range sampled here, 1900~7100 A, into five bins of equal 
bin-size in log The relationship of amplitude <r m versus rest-frame effective wavelength 
A is shown in Figured The mean value and scatter (standard deviation) of the sigma 
distribution in each wavelength bin is also overplotted. The result shows clearly a trend 
of increasing variability amplitude moving towards short wavelength for NLSl-type AGNs, 
similar to BLSl-type. It shows clearly that at a given wavelength NLSl-type AGNs have 
systematically lower variation amplitudes compared to BLSl-type, albeit large scatters. 
We test the significance of the differences in the a m distributions between the two samples 



4 The rest wavelength bin, in A, are 1900, 2500, 3300, 4200, 5500, 7100. 
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using the K-S test, yielding the chance probability for the same distribution as 1.2%, 9.7%, 
2.9% for the shortest three wavelength bins, respectively. 

The trend in Figured is in statistical sense only, and it would be interesting to examine 
whether it also holds for variability in individual objects. For each object, we estimate 
the maximum magnitude difference over the entire light curve in the blue (u) and red (i) 
bands, respectively, defined as Am = m t 2 — rn t i, where tl and t2 are the epochs of the two 
observations and t2>tl. The result is shown in the left panel of FigureE) the dotted line 
indicats the 1:1 relationship where the variations in the two bands are identical in both the 
sign and amplitude. Two conclusions can be inferred from the figure. First, Am in both the 
blue and red bands always keep the same signs, meaning that the two bands vary (fading or 
brightening) coordinately; That is to say, the optical/UV continuum level moving up and 
down systematically, rather than seesawing around. Second, for all but a few objects the 
variations in the blue are indeed larger than those in the red band. This suggests that, even 
in individual AGNs, the brightness variation tends to go hand-in-hand with variation of 
the continuum slope, with larger variability in blue than in red. This is the same for both 
NLS1- and BLSl-type AGNs. 

A direct consequence of such wavelength-dependent variability is that the optical/UV 
continuum becomes redder (bluer) when the overall continuum level becomes fainter 
(brighter). This can be demonstrated explicitly in the right panel of Figure[8l where the 
changes in the color (A(m u — m,j)) are plotted against the changes of magnitude. As can 
be seen, for all but a few objects, of both BLS1- and NLSl-type AGNs, the spectra redden 
(A(m u — rrii) > 0) while objects fade away (Am > 0), and vice verse. Furthermore, there is 
a strong correlation between the change of color (A(m u — mi)) and the change of magnitude 
(r s =0.79, P < 10~ 5 , the Spearman rank correlation test). This implies that the larger 
variation of the continuum level is, the more change of the spectral shape in general. We 
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note that, regarding the dependence of variability on wavelength and spectral shape, NLS1- 
and BLSl-type AGNs show very similar behavior. 



4.3. Time dependence of variability 



To characterize t he dependence of variati o n on time lags, structur e function is 



widely used (see, e.g., ICollier fc Peterson! l200ll ; IVanden Berk et al.ll2004l ). In this paper 



we construct the structure function with maximum-likelihood method, as described in 
Section[3l to parameterize the intrinsic variability amplitude for the ensembles instead of 
individuals. Firstly, for each object, we calculate the magnitude difference between any two 
observations separated by an interval, At, in the object rest frame, Am = m(t) — m(t — At). 
The error of Am is estimated using error propagation. The combined data points for all 
objects in the samples are then grouped into various bins according to time lags. The time 
lag bins are divided in a way that there are at least 600 data points in each bin. In each 
bin of time lag, the mean of the Am distribution is expected to be around zero for a large 
number of data points, since objects brighten or dim randomly and evenly in the statistical 
sense. The dispersion of Am distribution, in fact, reflects the degree of variations (at the 
given time lag). We take the intrinsic dispersion of the Am distribution in each bin as a 
measure of the 'averaged' amplitude of variability, i.e., the structure function. The intrinsic 
dispersion can be obtain with the assumption that Am distribution in each bin can be 
approximated by a Gaussian, which is verified in the majority of the time lag bins. The 
uncertainties of the intrinsic dispersion can also be obtained at the same time. In this way 
we construct the structure function with maximum-likelihood method. 

We first construct the structure function for NLSl-type and BLSl-type AGNs in 
each SDSS photometric band, for the sake of easy comparisons with previous results of 
BLAGNs and quasars. The derived structure functions are shown in Figure[9l in only g and 
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r bands for clarity. We find that the variability amplitude of NLSl-type AGNs, similar to 
BLSl-type, does increase with time lag from days to years in all of the SDSS bands. At 
a given timescale, NLS1- have systematically lower variation amplitudes than BLSl-type 
AGNs. A simple power law model is used to characterize the structure function at time 
lags greater than 10 days, 



where Atq and /3 are to be determined. The fitting results are given in Table[3] and Figure[9j 
It can be seen that a power law model can represent the structure function reasonably well. 
The inferred power law slopes for NLSl-type AGNs range from 0.32±0.02 to 0.40±0.01, 
mostly consistent with those of the BLSls within the uncertainties. 

Considering the redshift effect, a more physically meaningful approach is to construct 
the structure function in rest-frame wavelength bins instead of in photometric bands. We 
divide the whole rest wavelength range into five bins with equal bin-size in logarithm, as 
described in Section l4~2l The constructed structure functions for both NLSl-type AGNs 
and BLS1- in each bin are shown in Figure[TUJ As expected, the time lag dependence of 
variability presents in all the bins. Wavelength dependence variability is also clearly present 
in almost all of the time lag bins, which is consistent with what was found above and also 
in previous work for BLAGNs. Again, the two types of AGNs show remarkable similarity 
in the shape of their structure function on timescales longer than about 10 days. 

On shorter timescales (below about 10 days), the structure functions of the two types 
seem to differ, however. The structure function of NLSl-type AGNs continues to drop, 
whereas that of BLS1 seems to flatten towards short time lags, with some residual variability 
of several percent. This is the case for all the wavelength bins except for the shortest 
wavelength, for which the photometric uncertainties are the largest. Also we note that 
the wavelength- dependence effect becomes weak or even vanishes on such timescales. This 




(2) 
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may indicate a different mechanism of variability starts to become dominating. However, 
this is not the case for NLSl-type AGNs, of which the structure function continues to 
drop towards short time lags. Given the sparse sampling of the data used here, further 
confirmation is needed to verify this trend. 

The structure functions in various rest frame wavelength bins are also fitted with a 
power law model within the time lag range of greater than 10 days. The fitted curves are 
shown in Figured]] for the bin of 2500-3300 A as an example, and the results are given in 



Table[3j Alternatively, an exponential model (IMacLeod et al.l l2Q10h is also used 



SF(t) = SF^l - e -l^*lA)V2 ( 3 ) 

where SF^, and r, the time scales, are to be determined. With inferred SF^ we can 
estimate the long-term standard deviation of the variability as 0.5(S'-F oo ) 2 . It can be 
seen that, although both the functions can approximate the observed structure functions 
reasonably well (at r >10 days), we note that the exponential function gives a somewhat 
better description than the power law for the NLS1 sample, based on the xt statistic. 
This may be ascribed to a possible flattening of the structure function at the largest time 
lags, which appears to be more significant for the NLS1- than for BLSl-type AGNs (see 
FigurefTUI Figures fTT]) . However, the time lags in the current data set are not long enough 
to draw the conclusion. The inferred variability timescales are around one year and increase 
towards shorter rest-frame wavelength bins. The values of SF^, i.e., the maximum variation 
amplitude at the longest possible time lags, as expected, are larger at shorter wavelength 
than longer, and larger for BLS1- than NLSl-type AGNs in all the rest wavelength bins. 
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5. Discussion 

5.1. Implications to the optical/UV variability of NLSls 

The similarities found between NLS1- and BLSl-type AGNs in most of the properties 
of optical/UV variability, except for systematically smaller amplitudes of NLSls, indicate 
that their long-term variability must be driven by the same mechanism. Recently there 
are compelling evidences that the variability is intrinsic to AGN activity, specifically the 



detected correlations of optical/UV v ariabilities with BH masses and/or 



as investigated in our previous study 



flAi et al. 



Wilhite et al. 


2008; 


Bauer et al. 


2009; 


Zuo et al. 


2012) 



20101) and other research ( Wold et al. 



Iddington ratios. 



2007 : 



20121 ). The variability we observed may 



be caused by variations in the accretion rates [see Gaskell (2008) for a differe nt view. 



however], or some kind of (loc al) disturbance a nd its propagation in the disk (ICzerny 
20061 ). Assuming the first case, iLi fc Cad (120081 ) proposed a simple model of global change 
of the disk structure due to variation in accretion rates based on the standard disk, which 
explained the previously observed Mbh and luminosity dependence of variability. In fact, 
this model can also reproduce qualitatively an inverse amplitude-Lboi/^Edd relation (Li, 
private communication), which is MsH^dependent, though a quantitative comparison with 
our data is hampered due to the small sample size. 

Le t's consider the case of lo cal instability scenario. From the simple standard disk 



model ( Shakura fc Sunyaev 



19731 ) . the emission at a given radius R is dominated by a local 



blackbody with T e g peaking at a wavelength A, where 

hc/X ~ 2.8kT cS . 

The radiation is balanced by the energy generated per unit disk area at the radius R, 

aT c 4 ff ~ 3GM hh M/(87iR 3 ), 
where M is the mass accretion rate. Combining the two relations we find the radius R 
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where the local blackbody emission is peaked at A is, in units of the Schwarzschild radius 
R Sch , r = R/Rsch ~ (m/ M bh) 1/3 A 4/3 , where rh the mass accretion rate in units of the 
Eddington rate. This means that the scaled radius r, which dominates the emission in a 
given bandpass, enlarges with increasing Eddington ratio. If the disturbance is generated 
in the inner part of the disk and propagating outward, it would attenuate in the course of 
propagation. This would account for the ensemble weaker variability of NLSl-type AGNs 
compared to BLSl-type AGNs, since the radius r responsible for the emission in the same 
bandpass is shifted further out as the Eddington ratio increases. On the contrary, if the 
disturbance is generated in the outer disk and is propagating inward, the disturbance must 
be amplified instead. 

Reprocessing of the X-ray emissions, originating close to the central black hole and 
known with strong and rapid variabil ity, into the optical /UV ban ds also can cause the 



optical/UV variability in AGNs (e.g. 



Czerny 



2006 



Gaskcll 



20081 ). Yet the multi-bands 



variability correlation studies in individual AGNs indicate that the dominant mechanism 
for the long-te rm AGN optical/UV varia bility can not be reprocessing of variable X-ray 



emission (e.g. 



Arevalo et al. 



2008 



20091 1 . Our result about the statistically weaker 
variability of NLS1- than BLSl-type AGNs in optical/UV bands supports the argument, 
since NLSl-type AGN s often show strong X-ray variability on both short and long timescales 



e.g. 



Papadakis 



20041 ). However, on timescales of less than 10 days, it is likely that in this 



regime the X-ray reprocessing starts to take over and dominate the optical /UV variabili ty. 



as the relevant timescales involved in this process are generally very short (I Czerny 



2006|). 



We find that the variability component at short-term timescales seems to behave 
differently from the long-term one in that the wavelength-dependence is absent, or at least 
weak. And the structure function of BLSl-type AGNs begins to flatten at time lags less 
than 10 days. Interestingly, this flattening is not seen in the structure functions of NLSls, 
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which remain declining in the power law form. The effect of contaminations from emission 
line v ariability is excluded since the emission lines are m ostly not variable on timescales of 



days ( iVestergaard fc Peterson 



2005 



Wilhite et al 



20051 ). We propose that such a differe nce 



Uttlev 



may b e explained in terms of the X-ray reprocessing scenario. As suggested in 
( 120061 ). the interaction between optical/UV and X-ray variability maybe different between 
NLS1- and BLSl-type AGNs. The optical/UV emission regions for NLSl-type AGNs 
are outer in the accretion disk compared to BLS1-, which subtend only a smaller solid 
angle as seen from the X-ray sources; while in BLSl-type AGNs the optical/UV emissions 
originate much closer to the X-ray emitting region, which means much tighter interactions 



NLSl-tv 


pe 


Komossa 


2008 



AGNs appears to be generally weaker and softer than BLSl-type AGNs (see 
and references therein), which may also lead to generally weaker optical emission from the 
reprocessing process. Thus, the weakness or absence of this excess variation component 
may also contribute, at least partly, to the overall lower variability amplitudes in NLS1- 
compared to BLSl-type AGNs. 

If the dominant mechanism for the long-term optical/UV variability is related to 
accretion disk instability, we would expect that the characteristic variability timescales 
inferred from structure function are different between NLS1- and BLSl-type AGNs. The 
typical physical timescales associated with the accretion disk scale with the black hole mass. 
However, in our analysis the variability timescales inferred from the exponential fits of the 
structure function do not show statistically significant difference between the two types of 
AGNs given their large uncertainties. This may be ascribed, at least partly, to the relatively 
small range of black hole masses i n our sample. We note that in the analysis of SDSS 



Stripe-82 data for massive quasars, 



MacLeod et al. 



(120101 ) find the characteristic timescale 



correlates with the black hole mass with a power law index of 0.21 ±0.0 7. While AGNs with 
relatively smaller black hole mass, such as the NLSl-type AGNs are not included in their 
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sample. The property of the characteristic timescales of AGN variability may be examined 
by using large samples with a wider black hole mass range and monitored over longer 
time spans than the currently available data, such as the upcoming Large Synoptic Survey 
Telescope. 



5.2. Implication for the nature of NLSl-type AGNs 

We have demonstrated that, on all timescales observed so far, NLSl-type AGNs 
show systematically less optical/UV variations compared to BLSl-type AGNs. This is in 
agreement with the de pendence of th e variability on the Eddington ratio for BLAGNs found 



in our previous paper (lAi et al. 



20101 ). In fact, as we argued, the optical/UV variability may 
be considered as a new EV1 parameter, with NLSls lying at one extreme end. This property 
gives interesting insight into the nature of NLSl-type AGNs. Currently there are three 
scenarios regarding the nature of NLSl-type AGNs. The widely acce pted model invokes les s 
roles accreting at high rate s compared to BLSls (e.g. 



massive black 



Sulentic et al. 



2000 



Wang k, Netzer 



Mineshige et al. 



2000 



20031 ). Alternatively, the narrow Balmer line width 



can also be explained by a disk-like low-ionization BLR close to face-on (the "orientation 
scenario", e.g., Osterbrock & Pogge 1985, Collin & Kawaguchi 2004), or by a more distant 
BLR compared to normal BLSls with a similar black hole mass or nuclear luminosity 
("distant BLR scenario", Wandel & Boiler 1998). Although face-on inclination is found in 
a small number of extreme radio-loud NLSls showing blazar-like properties (Yuan et al. 
2008, abdo et al. 2009, Gu & Chen 2010), comparative studies of radio properties between 
NLS1 and BLS1 disfavor the orientation scenario in general (Zhou & Wang 2002; Komossa 
et al. 2006; Zhou et al. 2006). Other lines of evidence, such as polarization properties 
and correlations involving narrow-line luminosity, suggest that orientation at most plays a 
secondary role in explaining NLS1 phenomenon (Komossa 2008 and references therein). 
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Neither the "orientation" nor "distant BLR" scenario can explain naturally the 
extension of EV1 correlations involving variability. For the "distant BLR" scenario, it is 
hard to link the size of the BLR to the optical/UV continuum variability, that is local 
to the accretion disk. Also, since the optical continuum emission is believed to come 
from an accretion disk, as long as the disk is directly seen, variability amplitude should 
not be dependent on the orientation of the disk. Optical/UV variability may be caused 
by obscuration by the outer part of the disk; however, a lack of periodic variability 
corresponding to the disk rotation has ruled out this possibility, (e.g., Peterson & Bentz 
2006; Wold et al. 2007). In comparison, the commonly accepted model for NLSls with 
smaller Mbh and higher Lboi/-^Edd is favored. Their weak variability fits more naturally 
into the scenario that variability amplitude is governed by Z/boi/^Edd, as suggested in our 
previous paper. 

6. Summary 

We present the ensemble variability property of a sample of NLSl-type AGNs that 
have multi-epoch photometric observations in the SDSS Stripe-82 region. For NLSl-type 
AGNs the variability was hitherto poorly explored as a class due to both the limited sample 
sizes and few data available. As a direct comparison a control sample of BLSl-type AGNs 
is also compiled. We introduce a novel parametric parametric maximum-likelihood method 
in variability studies, to quantify intrinsic variability amplitudes and, for the first time, 
their confidence intervals, which improves upon the commonly adopted methods in the 
literature. This method makes full use of the measurement error of each data point, and 
has the advantage of being sensitive to low-level variations. 

We find that the majority of the NLSl-type AGNs are in fact variable on timescales 
from about several days to a few years, e.g., more than 80% of the objects varied at 
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levels > 0.05 mag (standard deviation) in the SDSS u and g bands. NLSl-type AGNs 
have systematically smaller variability compared to BLSl-type AGNs. This is consistent 



with the previously 



flWilhite et al. 



2008 



bund anti-correlation of AG^ 



Ai et al. 



2010 



MacLeod et al. 



varia bility with the Eddington ratio 



20101 ) . We present the light curves and 



some of the key parameters for 22 NLSl-type AGNs that show relatively large variations 
(standard deviation a > 0.1 mag in the r-band), which may deserve further monitoring 
observations. It is also found that, similar to BLSl-type AGNs (including quasars), 
the optical/UV variability is wavelength-dependent — the shorter wavelength, the larger 
variations. We derive, for the first time, the ensemble structure function of variability for 
NLSl-type AGNs that span nearly three orders of magnitudes in time lags (from a few days 
to ~ 10 3 days). The structure functions of NLS1- and BLSl-type AGNs show remarkably 
similar profiles on time lags > 10 days, with the former having smaller amplitude at a 
given time lag. The structure function can be approximated by either a power law or an 
exponential function, though the latter seems to give a better description for NLSl-type 
AGNs. On time scales of < 10 days, the variability amplitudes of BLSl-type AGNs remain 
somewhat higher than those extrapolated from the power law structure function, indicating 
the presence of some excess variations. However, this is not the case for NLSls. We suggest 
this excess variation component to be likely arising from reprocessing of the X-rays, which 
might be more important in BLS1- than in NLSl-type AGNs. 

In conclusion, NLS1- and BLSl-type AGNs share very similar properties of long-term 
optical/UV variability in many ways, such as the wavelength and time lag dependence of 
variability. This suggests that most likely the same physical mechanisms are at work to 
produce the long-term optical/UV variability in both types, in spite of the possibly different 
modes of their accretion flows. Future photometric monitoring programs for larger samples 
of AGN spanning wide parameter ranges, with a large dynamic range of sampling rate, are 
needed to address the questions related to the optical/UV variability, such as the behavior 
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at time lags less than a week and the characteristic variability timescales. 
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Fig. 1. — Redshift and i-band luminosity distribution of our NLSl-type AGN (filled circles) 
and BLSl-type (open circles) samples. The dashed curve represents the limiting magnitude 
for optically-selected quasars in the SDSS (m^ = 19.1). 
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Fig. 2. — Distributions of the H/3 linewidth, relative strength of Fe II emission, black hole 
mass and Eddington ratio for the NSL1- (solid) and BLSl-type (dash-dotted) AGNs. 
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Fig. 3. — Example of the best estimates and confidence contours of the mean (m) and 
dispersion <x m of the intrinsic magnitude distribution at the 90% and 60% confidence levels 
in ugriz bands for one NLSl-type AGN. 



-35 - 



0.25 

0.2 o - 

^ 0.15 - 

en 
o 

E 

e 
t> 

0.10 - 
0.05 - 




0.00 I I , , , , ■ i i i i , , , , ■ i i i i 

1 10 100 

xl 



Fig. 4. — The r-band variability amplitudes versus reduced chi-squares for NLSl-type AGNs. 
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Fig. 5. — Accumulative distribution of the variability amplitudes, a m , for NLS1- (filled 
circles) and BLSl-type AGNs (open circles) in the five SDSS photometric bands. The x-axis 
ranges truncate at 0.3 mag for clarity. 
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Fig. 6. — Light curves of the NLSl-type AGNs in r band with a m > 0.1 mag. 
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FigureE] — continued 
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Fig. 7. — Variability amplitude as a function of rest-frame wavelength for NLS1- (left) and 
BLSl-type AGNs (right). The mean value in each bin is also shown with the error bar 
representing the standard deviation of the distribution in each bin. 
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Fig. 8. — left: Comparison of the maximum magnitude differences between the u and % bands 
for each NLS1- (filled circles) and BLSl-type AGNs (open circles); right: the changes in the 
color versus the maximum magnitude differences. 
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Fig. 9. — Structure functions of the NLS1- (filled circles) and BLSl-type AGNs (open circles) 
in g (blue) and r (red) bands. The errors are at 90% confidence level. The dashed lines are 
the power law model fittings to data at At >10 days. 
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Fig. 10. — Structure functions of the BLS1- (left ) and NLSl-type AGNs (right) in rest-frame 
wavelength bins of AA(1900-2500 A) (blue), AA(2500 -3300 A) (green), AA(3300 -4200 A) 
(red), AA(4200 - 5500 A) (pink), and AA(5500 - 7100 A) (black). 



-43- 



i i i i n-i-f 




10 100 1000 

Rest Frame Time Lag, At (days) 



Fig. 11. — Structure functions of the NLSl-type AGNs (filled circles) and BLSl-type AGNs 
(open circles) in rest wavelength bin, 2500-3300 A. The fittings with power law (dotted) and 
exponential model(dashed) are also presented. 
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Table 1: Fraction of variable objects and averaged magnitude of variations. 



NLSl-type AGNs BLSl-type AGNs 



bands P(x 2 \ v) < 0.001 a 


a m > 0.05 6 


a m > 0.10 c 




P(X 2 1 v) < 0.001 Q a m > 0.05 b o 


m > 0.10 c 




u 85.4 


89.1 


53.4 


0.11 


93.5 90.6 


71.5 


0.14 


g 98.3 


85.4 


44.8 


0.10 


99.1 95.3 


60.1 


0.13 


r 98.3 


76.2 


41.3 


0.09 


99.1 90.6 


56.9 


0.12 


i 100. 


74.5 


39.6 


0.08 


97.2 89.7 


49.6 


0.10 


z 81.3 


69.4 


22.4 


0.07 


85.0 80.3 


38.2 


0.09 


a Fraction of sources with P(x 2 


| u) < 0.001; 


''Fraction of 


sources with 


&m > 0.05; c Fraction of sources with 


a m > 0.10; 





d Mcan a m in each band 
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Table 2: NLSl-type AGNs with a r > 0.1 mag 



SDSS Name 


Redshift 


FWHM(H/3) 


lo£?"A Tj \ c; i nn 
* A 5 1 U U 




R4570 


(r) 




' 9 










mavA g 


may A r 


(1) 


(2) 

V-J 


(3) 






(5) 


(6) 




(7) 










(9) 


do) 


T003431 74-001 SI 2 7 


0.381 


1701+31 


44.47 





27+0.01 


18.14 





11+0 


05 





10+0 


04 


0.37 


0.33 


1010737 01-00191 1 6 


0.737 


1853+57 


45.10 





78+0.04 


18.34 





19+0 


08 


o 


18+0 


,08 


0.56 


0.59 


J011712. 82- 005817.4 


0.485 


1 987+65 


44.51 





49+0 03 


19.24 


o 


15+0 


05 


o 


14+0 


04 


0.49 


0.47 


T01 2824 21 4-001 925 2 


0.419 


1 K r=o+ 1 32 


44.13 







19.70 


o 


16+0 


05 


o 


17+0 


05 


0.55 


0.53 


T01 3509 51 4-000252 


0.744 


1890+129 

i u -"7 v_y i Lie 


44.73 





96+0 12 


18.99 


o 


17+0 


07 


o 


15+0 


06 


0.49 


0.43 


10,9^0,1 1 8-1-00,1 7A^ A 


0.359 


i on 


44.04 


X 


07+0 14 


19.41 


o 


17+0 


05 


o 


11+0 


03 


0.59 


0.40 


102^910 3Q-00223Q 8 


0.360 


208^+^2 


44.52 







18.20 


o 


16+0 


06 


o 


15+0 


06 


0.51 


0.40 


T030335 764-004145 


0.669 


i qqq-|-a'3 

J- ■ ' ■ ' 1 UP 


44.87 





i q-l-n (12 


19.24 


o 


15+0 


04 


o 


14+0 


01 


0.59 


0.57 


J031022. 114-004130.0 


0.655 


1972+203 


44.60 


x, 


44+0.23 


19.74 


o 


12+0 


03 


o 


10+0 


03 


0.46 


0.40 


J031427.46-011152.3 


0.387 


1709+22 


44.59 


0. 


50+0.01 


18.15 


o 


14+0 


05 


o 


14+0 


06 


0.54 


0.50 


J031532.28+005503.6 


0.487 


2059+165 


44.12 





00+0.05 


19.74 





,23+0 


08 





,24+0, 


09 


0.98 


0.87 


J031630.78-010303.6 


0.368 


1174+87 


44.28 





79+0.07 


18.73 





.15+0 


04 





,12+0 


03 


0.55 


0.41 


J034233.05-001448.5 


0.695 


1502+118 


44.54 


0. 


53+0.18 


19.52 





,16+0 


06 





17+0 


06 


0.49 


0.54 


J204527.70-003236.2 


0.296 


1473+53 


44.20 


1. 


13+0.04 


18.37 





,14+0 


06 





,11+0 


05 


0.44 


0.38 


J204910.97+001557.3 


0.362 


1556+41 


44.17 





25+0.03 


18.80 





,12+0 


06 





10+0 


01 


0.47 


0.42 


J212130.14-000527.1 


0.585 


1898+39 


45.32 





95+0.03 


17.39 





,15+0 


07 





,13+0 


05 


0.40 


0.44 


J212341. 16+010557.9 


0.389 


1473+52 


44.23 


0. 


75+0.07 


19.33 





,16+0 


06 





14+0 


05 


0.66 


0.57 


J214743.56+004028.2 


0.643 


1987+72 


44.94 





80+0.03 


18.98 





,22+0, 


08 





,22+0, 


08 


0.68 


0.61 


J221451. 27+000725.1 


0.573 


1780+87 


44.69 


1 


01+0.08 


18.93 





,09+0 


01 





11+0 


05 


0.32 


0.37 


J222732.93+010912.3 


0.546 


1472+46 


44.61 





68+0.04 


19.11 





,16+0, 


06 





,15+0 


05 


0.58 


0.53 


J230555.49+005946.2 


0.719 


2161+47 


45.17 





44+0.02 


18.14 





,11+0 


04 





,11+0 


05 


0.33 


0.35 


J233517.62+010307.7 


0.624 


2095+109 


44.66 





01+0.03 


18.84 





,10+0 


03 





10+0 


03 


0.51 


0.50 



Note. — Col.(l): SDSS name; Col. (2): redshift; Col. (3): H/3 lincwidtli (kms^ 1 ); Col.(4): monochro- 
matic luminosity at 5100 A (ergss -1 ); Col. (5): intensity ratio of the Fe II multiplets to H/3; Col. (6): mean 
magnitude in r band; Col. (7)- (8): best estimated variability amplitudes and corresponding errors at 90% 
confidence level from maximum-likelihood method; Col.(9)-(10): maximum magnitude changes 
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Table 3: Results of model fits to the structure functions at At > 10 days 

Band NLSl-typc AGNs BLSl-typc AGNs 



Power Law 









Aro(d) 


P 


Aro(d) 





u 






7.3xl0 4 


0.33±0.01 


1.7x10* 


0.36±0.01 


g 






6.3xl0 4 


0.36±0.01 


2.2x10* 


0.37±0.01 


r 






5.4x10* 


0.40±0.01 


3.1x10* 


0.35±0.01 


i 






8.1x10* 


0.39±0.01 


4.3x10* 


0.36±0.01 


z 






3.6xl0 5 


0.32±0.02 


5.7x10* 


0.37±0.01 


AA 1900- 


■2500 A 


1.4xl0 6 


0.30±0.02 


2.1x10* 


0.34±0.01 


AA25OO- 


3300 


A 


6.6x10* 


0.36±0.01 


2.9x10* 


0.35±0.01 


AA330O- 


4200 


A 


5.7x10* 


0.38±0.01 


2.6x10* 


0.36±0.01 


AA420O- 


5500 


A 


6.5x10* 


0.40±0.01 


3.2x10* 


0.38±0.01 


AA550O- 


7100 


A 


3.7xl0 5 


0.33±0.02 


8.6x10* 


0.34±0.01 


Exponential 








r(d) 




r(d) 


SFoo 


AA1900- 


2500 


A 


217±28 


0.21±0.01 


282±24 


0.32±0.01 


AA2500- 


3300 


A 


336±33 


0.21±0.01 


309±21 


0.28±0.01 


AA3300- 


4200 


A 


404±40 


0.21±0.01 


343±21 


0.28±0.01 


AA4200- 


5500 


A 


562±58 


0.20±0.01 


424±27 


0.26±0.01 


AA5500- 


7100 


A 


410±67 


0.14±0.01 


401±39 


0.21±0.01 
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A. Re-calibration of Supernova Survey data 

Re- calibration of SN data is performed in five bands separately in a field-by-field 
(lOOarcmin 2 patches) manner, and only fields with "good" image quality (flagged as "good" 
or "acceptable") were used for photometric calibration or light curve construction in our 
work. Photometric conditions can be regarded to be unchanged across the field. For each of 
the SDSS SN survey fields, we only use the most overlaying DR5 field to do calibration. In 
the case of more than one overlapped DR5 fields, we use the one with the larger or largest 
overlapping area. In each field and each band, 'standard stars' with good measurements 
(high quality photometry) are selected with the following criteria: 

1. Sources must be 'star'. 

2. The processing flags BRIGHT, SATURATED BLENDED, EDGE, NOPROFILE, 
INTERP .PROBLEMS, and DEBLEND_PROBLEMS are not set in any bancfl. 

3. Sources must be brighter than the magnitude limits which are 22.0, 22.2, 22.2, 21.3, 
20.5 in ugriz bands. 

4. Potentially variable sources were excluded, including QSOs, kno wn variable star s, 



etc. Here QSO identification was based on SDSS quasar catalog. III. (ISchneider et al. 



200 5]), and varia 



jle sta rs based on the General Catalog of Variable Stars, 4th Ed. (GCVS4) 



(IKholopov et al. 



19981). 



The selected stars were then cross-identified between the two most overlaying fields 
with a matching radius of l". Any stars with no match or more than one matches within 
this radius were discarded. We computed various photometric statistics, such as mean 
magnitudes and mean standard errors (weighted by statistical errors), and rms scatter of 



Sec http:/ / www.sdss.org/dr7/products/catalogs/flags.litml 
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the magnitude difference of the matched stars. Stars with magnitude difference more than 
two times of the rms scatter were rejected in order to avoid variable stars, mismatched 
stars, or bad data points, etc. The above procedure was repeated until no more star was 
excluded. The final number of calibration stars is typically ~ 50 — 200 for each pair of fields, 
which is large enough to do calibration with an acceptable uncertainty. For each of the 
fields, the differential magnitudes of the calibrating stars between the SN survey and DR5 
observations (Am = msN — ^drs) are calculated, and their weighted mean is set to zero, 
(Am) = 0. In this way the zero-point offsets of the SN survey photometry are determined, 

The internal uncertainties of the calibrated zero-points of the SN survey photometry 
with respect to the reference DR5 data, as measured by the standard deviation of the 
weighted mean of Am, i.e., <7(Am}> are less than 0.01 mag for all the fields. However, the 
SDSS SN survey data are calibrated using different fields of SDSS-I DR5 data, which 
have typical systematic photometric zero-point errors (<x^,f 5 ) of 0.01 mag for the g, r, and 



% bands, 0.02 mag for the z band, and 0.03 mag in the u band ( llvezic et al 



20041 ). Thus 



the total systematic uncertainties of the SN survey photometry are estimated to be the 



combination in quadrature of the above two terms, i.e. J (e>"(Am}) 2 + (o^f 5 ) 2 - The overall 
photometric errors are then a contribution of both the systematic and statistical errors. 
The latter are given in the SN survey catalogs and range from a few percent to 0.1 or 
higher toward faint magnitudes, depending on the bands. The final photometric errors of 
the calibrated SN survey data have a median of ~ 0.03 mag in the g, r, and i bands, and 
somewhat higher (« 0.04 mag) in the u and z bands. They are comparable to those of the 
DR5 data observed in SDSS-I. 

As an example, the SN data of one field before and after calibration are illustrated in 
Fig. [T2J Actually our photometric accuracy is much better than the rms scatter shown in 
Fig. [121 because most of the objects in our samples are brighter than 20.0 mag in r band. 
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Fig. 12. — Illustration of the zero-point calibration for one field in SDSS Supernova survey, 
with the cross-identified star from the most overlaying filed in DR5. The dots in the left 
panels represent magnitude differences between the DR5 and the original SN magnitudes, as 
a function of the DR5 magnitude. The dashed line shows the mean magnitude differences 
weighted by statistical errors. The middle panels show the magnitude differences after our 
calibration. The histograms of the magnitude differences before (dashed line) and after (solid 
line) calibration are compared in the right panels. 
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B. Amplitude of Variability and the Maximum- likelihood method 



A common practice as used in the literature is to estimate the contribution of the 
errors of measurement to the observed sca tter, and then subtract it from the latter 



( jVanden Berk et al 



2001 



Sesar et al. 



20071 ). The observed rms scatter E is calculated as 



1 N 



(Bl) 



where rrii are observed magnitudes of N observations and (m) is their weighted mean. The 
intrinsic variation amplitude a m can be estimated as 



^method— 1 



0, otherwise. 



(B2) 



where £ is the term representing the amount of scatter caused by measur emental errors. T 



no 



Sesar et al. 



fl2007h 



estimation of £ has some variances among the literature. For example 
used the mean photometric errors of SDSS observations as a function of magnitude which 
was fitted from the data assuming most stars are not variable; a theor etically expected 



Strateva et al. 



(j200lh . Here we estimate 



£(m) as a function of magnitude was also given in 
£ as the mean square value of the errors £j associated with the individual magnitude rrii, i.e. 

1 N 

-Ye 



(B3) 



1=1 



This estimation is similar that used in Rodiguez-Pascual et al. (1997), which was in the 
flux domain rather than magnitude though. 

In this paper, we introduce a new method — the maximum-likelihood method — to 
quantify a m and its uncertainty, in the line of the parametric statistics approach. There are 
two assumptions involved in this method. Firstly, we assume that, for a given measurement, 
the probability function of the random fluctuations in the measurement follows a Gaussian 



-51 - 



distributiorjj, with a standard deviation of a r , where a r is the photometric error for this 
measurement. Furthermore, for simplicity, we approximate the intrinsic distribution of the 
magnitude of a variable object as a Gaussian^, with a mean (m) and a standard deviation 
a m - both parameters are to be estimated. cr m gives a measure of the amplitude of variability, 
and is expected to be consistent with zero in case of no variation. 

The final probability function of observing an object to have a magnitude m with a 
photometric uncertainty is then 



p(m) = exp 



[m — (m))' 



(B4) 



2<p 2 

Here cp = {cr^ nethod _ 2 + of) 1 ^ 2 - Suppose, for a given object, a set of values for m are 
observed, mi, m^, with the corresponding (different) accuracies of the measurements 
which are characterized by the standard deviations a r i, a r N- The likelihood function, 
L((m), (Trnethod-2), of this data set is 



N 

TT 1 (m,i — {m)) 2 

A 1 +(7 2 U/2 eX P 



2 ( <T meihod-2 +<J ri) 



The maximum-likelihood estimates of (m) and cr me thod-2 for a given object are obtained by 
minimizing S = — 21nL 

This method can also be used to find the confidence intervals for the interesting 
parameters, by applying the standard A% 2 techniques to the above S-function (see 



6 This is supported by the fact that for stars (mostly non-variable) repeatedly observed 
in the SDSS with similar photometric errors, their magnitude distributions can be well 
described Gaussian (Ivezic et al. (2003)). 

7 We consider this as a reasonable approximation, as it is found from our data that, for 
some objects with a much larger amplitude of intrinsic variations than their typical photo- 
metric errors (i.e. the latter is negligible), their magnitude distributions can be reasonably 
well described by a Gaussian. 
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e.g. Avni 1976) . Suppose the parameter vector 9 has two components, and the first 
"interesting component" 6*j consists of q parameters ...,9 q ), which need to be es timat ed 
simul taneously without considering the second "uninteresting component" 9 U (see 



Avnil 



19761 ). The confidence interval at a probability level a for 9i can be computed by finding 



the set of all values of 9i such that 

S (9 i, minimize over 9 U ) — S min < A(q,a), (B5) 

where A(q,a) follow a \ 2 distribution with a degree of freedom q. In our case, there are 
two interesting parameters (q = 2), and the corresponding A(g,a)=2.30, 4.61, 9.21 for 
probability levels a = 68%, 90%, 99%, respectively. 

Intrinsic variability amplitudes estimated via the above two methods are compared in 
Fig. [13] and are found to be in excellent agreement. 
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Fig. 13. — Comparison of the estimated variability amplitudes with the two methods in r 
band. The other four bands give similar results. 



